The generalized non-linear Schrodinger model 

on the interval 



Anast DoikouJJ Davide Fioravantrr and Francesco Ravaninirr 

"66 
O 
o , 

. " University of Bologna, Physics Department, INFN Section 

■ Via Irnerio 46, Bologna 40126, Italy 



(D 
£^ 

m 
> . 

^ . Abstract 

in 

The generalized (H-l)-D non-linear Schrodinger (NLS) theory with particular inte- 
^3 I grable boundary conditions is considered. More precisely, two distinct types of bound- 

I ary conditions, known as soliton preserving (SP) and soliton non-preserving (SNP), are 

' implemented into the classical gl^ NLS model. Based on this choice of boundaries the 

relevant conserved quantities are computed and the corresponding equations of mo- 
^ I tion are derived. A suitable quantum lattice version of the boundary generalized NLS 

5^ I model is also investigated. The first non-trivial local integral of motion is explicitly 

computed, and the spectrum and Bethe Ansatz equations are derived for the soliton 
non-preserving boundary conditions. 



^ doikou@bo.infn.it 
^fioravanti@bo. infn.it 
^ravanini@bo. infn.it 



1 Introduction 



After various investigations (cf. below for detailed references) it is now well established that 
any integrable system (with finite or infinite degrees of freedom) based on the higher rank 
algebras gl^ or Uq{glN) may be endowed with two distinct types of integrable boundary 
conditions. These boundary conditions are known as soliton preserving (SP), traditionally 
recognized in the framework of integrable quantum spin chains (finite number of degrees 
of freedom) [l]-[8], and soliton non-preserving (SNP) originally introduced in the context 
of classical integrable field theories (infinite number of degrees of freedom) [9], and further 
investigated in [lOl [11], [12] . SNP boundary conditions have been also introduced and studied 
for integrable quantum lattice systems [6], [I3]-[T6], and their quantum integrability was 
first shown in [13]. From the physical point of view the SP boundary conditions oblige a 
particle-like excitation to reflect to itself: no multiplet changing occurs. The SNP boundary 
conditions, on the other hand, force an excitation to reflect to its 'conjugate', namely to 
an excitation carrying the conjugate representation. From the algebraic perspective the two 
types of boundary conditions are associated with two distinct algebras, i.e. the reflection 
algebra [H [2] and the twisted Yangian respectively pTl [T8] . 

The study of the underlying algebraic structures defined by the Yang-Baxter and reflec- 
tion equations is in general of great consequence, both at classical and quantum level, not 
only for integrable systems per se, but also for other relevant problems in theoretical physics. 
For instance, in the context of the AdS/CFT correspondence [I9] it is known that from the 
string theory side the relevant classical integrable model is a sigma model (see e.g. [201 El]), 
that is a field theory with infinite degrees of freedom. From the quantum gauge theory side 
on the other hand the loop contributions are apparently given by integrable quantum 1-D 
lattice models with a finite number of degrees of freedom [221 [23]. As a consequence, a cru- 
cial point would be the formulation of a discrete-quantum counterpart of the aforementioned 
classical sigma model. In this respect, the knowledge of the discrete-quantum Lax opera- 
tor would facilitate the derivation of the relevant Hamiltonian, and of the other charges 
in involution, as well as of the exact nested Bethe Ansatz equations. In fact, up to date 
only the asymptotic forms of the would-be- exact Bethe Ansatz equations are known (see e.g. 
[23]). Thus a rigorous derivation of these equations would undoubtedly be of great physical 
significance, as proven when corrections to the asymptotic regime are available [211 [25] . 

In the frame of classical continuum theories the SNP boundary conditions have been 
primarily investigated up to now [9]-[T2]. Therefore it is of great importance to further 
analyze the other set of boundary conditions, i.e. the SP ones within this context. In the 
present study we examine both SP and SNP boundary conditions for the classical generalized 

^They are defined as those boundary conditions that preserve the integrabihty of the system. 
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NLS model, and by using Hamiltonian methods [26] we derive the relevant integrals of 
motion, and also specify the corresponding classical equations of motion. Note that in [27] 
the quantum glj^ NLS model on the half line was studied, based on the reflection algebra 
(i.e. SP boundaries), primarily from the point of view of the underlying symmetry algebras. 
It should be stressed that although in most classical continuum theories the SNP boundary 
conditions have been analyzed (see e.g. [ni[IIl[I2|) this is the first time that such boundary 
conditions are implemented within the generalized NLS model. Here we consider the classical 
NLS model as a simple example, however our main motivation is to search for all possible 
boundary conditions in other classical theories such as affine Toda field theories, principal 
chiral models and others. From the viewpoint of quantum lattice models on the other hand 
the extensively analyzed boundary conditions are the SP ones, thus we focus here mostly 
on the SNP case for a lattice version of NLS. In particular, we consider a suitable lattice 
version of the NLS model [281 129] with SNP boundary conditions, and we derive the exact 
spectrum and the corresponding Bethe Ansatz equations. 

The outline of this article is as follows: in the next section we present the generic algebraic 
setting for classical models on the full line and on the interval. More precisely we introduce 
the classical Yang-Baxter equation and the underlying algebra for the system on the full 
line. In the case of a system on the interval we distinguish two types of boundary conditions 
based on the classical versions of reflection algebra (SP) and twisted Yangian (SNP). Next 
the NLS model on the full line is reviewed and an explicit derivation of the local integrals of 
motion by solving the auxiliary linear problem [26] is presented. In section 3 being guided 
by the same logic and adopting Sklyanin's formulation [2] we derive the integrals of motion 
of the gljsf NLS system on the interval with SP and SNP boundary conditions. Moreover, the 
corresponding classical equations of motion are obtained for various boundary conditions. 
In addition the usual NLS model with a reflecting impurity is investigated in the same 
spirit. In section 4 a suitable lattice version of the NLS model is investigated. After a brief 
review on the model with periodic boundary conditions we deal with the model with open 
boundaries. First the spectrum and Bethe Ansatz equations are derived for the usual lattice 
NLS. Finally, the SNP boundary conditions are considered for the generalized NLS system. 
The first non-trivial local integral of motion is explicitly specified for particular choice of 
boundary conditions, and the spectrum and Bethe Ansatz equations are deduced for the 
simplest boundaries. 
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2 The general setting 



The line of attack which will be adopted for the study of the gl^ NLS model with integrable 
boundaries is based on the solution of the so called auxiliary linear problem [26]. It is 
therefore necessary to recall at least the basics regarding this formulation. Let \1/ be a 
solution of the following set of equations 



— = U(x,t,A)vl/ 



'dt 



V(x,t,A)^ 



(2.1) 
(2.2) 



with U, V being in general n x n matrices with entries functions of complex valued fields, 
their derivatives, and the spectral parameter A. The monodromy matrix from ( 12. ip may be 
then written as: 



T{x,y,X) = Vexp\ / lS{x' ,t, X)dx' 



(2.3) 



The fact that T is a solution of equation (12. ip will be extensively used subsequently for 
obtaining the relevant integrals of motion. Compatibility conditions of the two differential 
equation (12. ip . (12. 2p lead to the zero curvature condition 



U- V' + 



U, V 



0, 



(2.4) 



giving rise to the corresponding classical equations of motion of the system under consider- 
ation. 

There exists an alternative description, known as the r matrix approach (Hamiltonian 
formulation). In this picture the underlying classical algebra is manifest in analogy to the 
quantum case as will become quite transparent later. Let us first recall this method for a 
general classical integrable system on the full line. The existence of the classical r-matrix, 
satisfying the classical Yang-Baxter equation 



ri2(Ai - A2), ri3(Ai) + r23(A2) + ri3(Ai), r23(A2) 



0, 



(2.5) 



guarantees the integrability of the classical system. Indeed, consider the operator T{x,y, A) 
satisfying 

|Ti(a;,y,t, Ai), r2(a;,2/,t, A2)} = ri2{\i - X2), Ti{x,y,t, Xi)T2{x,y,t, X2) ■ (2.6) 
Making use of the latter equation one may readily show for a system in full line: 

lntr{T(x,2/,Ai)}, lntr{T{x,y, X2)}] = (2.7) 
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i.e. the system is integrable, and the charges in involution -local integrals of motion- may 
be attained by the expansion of the object lntr{T{x,y, X)}, based essentially on the fact 
that T also satisfies (12.11) . 

Our main aim here is to consider the gl^ NLS model on the interval. For this purpose 
we follow the line of action described in [26], but using Sklyanin's formulation for the system 
on the interval or on the half line (see also [30] for the sine Gordon on the half line). We 
briefly describe this process below for any classical integrable system on the interval. In this 
case one has to derive a modified transition matrix T, based on Sklyanin's formulation and 
satisfying the following Poisson bracket algebras M, T, i.e. classical versions of the refiection 
algebra and twisted Yangian respectively. It will be convenient for our purposes here to 
introduce some useful notation. Let 

r*i2W = ri2(A) for SP, r^X) = fi2(A) = Vi r'^^h^) Vi for SNP 

r(A) = T~\-X) for SP, r(A) = V T*(-A) V for SNP 

"1/ = antid(l, !,...,!) or V = antid{i, — z, . . . , — z) for n even only (2.8) 

In general V can be any matrix such that V"^ = I, for instance V = Iis also a possible choice 
(see e.g. [9]). Then the defining relations describing the classical refiection algebra and the 
twisted Yangian respectively, may be written in the following compact fornj^: 

|Ti(Ai), r2(A2)}= [ri2(Ai-A2), ri(Ai)T2(A2)" 

+Ti(Ai)rt2(Ai + A2)r2(A2) - r2(A2)r*2(Ai + X2)r,{X,). (2.9) 

To construct the generating function of the integrals of motion one also needs c-number 
representations of the algebra M (T) satisfying (12. 9p for SP and SNP respectively, and also 



Kf(Ai), A'2^(A2))=0. (2.10) 

The modified transition matrices, compatible with the corresponding algebras M, T are given 
by the following expressions [2]: 

r{x,y,t,X)=T{x,y,t,X) K~{X) f{x,y,t,X) (2.11) 

and the generating function of the involutive quantities is defined as 

t(x, y, t, X) = tr{K+{X) T{x, y, t, A)} (2.12) 



■'Note that the classical versions of the reflection equation and the twisted Yangian arc provided in general 
by more involved expressions for generic r matrices. In the present study we focus on r matrices satisfying 
''12(A) = r2i(A) (fi2(A) = r2i(A)), and in this case (|2.9p . are valid. 
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Due to ([23D it can be shown that flXTOj) 

\^t{x, y, t, Ai), t{x, y, t, A2)} =0, Ai, A2 G C. (2.13) 

Technical details on the proof of classical integrability are provided in Appendix A. 

By expanding In t{X) in powers of A~^ one recovers the local integrals of motion of the 
considered system, and this is achieved in the subsequent sections. Among the local integrals 
of motion there exist naturally the Hamiltonian, which also provides information regarding 
the corresponding equations of motion. This is in fact the formulation we are going to 
assume for the NLS system on the half line, although an alternative strategy would be to 
derive the modified Lax pair, compatible with the boundary conditions chosen, and hence 
the associated equations of motion (see e.g. [9]). Nonetheless, the rigorous derivation of 
the modified Lax pair is essentially based on the existence of local integrals of motion |26] , 
therefore the viewpoint adopted here is arguably the most natural. 



3 Classical local integrals of motion 

The main objective in this section is to solve the auxiliary linear problem for the generalized 
NLS model on the interval. Before however we proceed to the study of the model on the 
interval we shall briefly review the system on the full line. In any case, these results will be 
relevant for the boundary well. 



3.1 The generalized NLS on the full line 

We shall hereafter focus on the gl^ NLS model. Consider the classical r matrix 

P ^ 
r(A) = - where P = V Eij ® Eji (3.1) 
A ' ^ 

P is the permutation operator, and {Eij)ki = SikSji. The Lax pair for the generalized NLS 
model is given by the following expressions [26] : 

U = Uo + AUi, V = Vo + AVi + A^Va (3.2) 

where (see also [3T] ) 

N~l N-l 

Ui = — (^ Eii - Enn), Uo = y/K^^i^piEiN + tpiEm) 

i=l i=l 
N-l N-l 

Vo = m ^ {^JiipjEij - \^i\^ENN) - iy/K^i^-EiN - ipiEm), 

i, j=l i=l 

Vi = -Uo, V2 = -Ui (3.3) 
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and ipi, ipj satisfja: 

M^)^ i^j{y) \ = \ '^li.x)^ i^jiy) = O, \ tpi{x), ipj{y) \ = Sij 6{x - y). (3.5) 



From the zero curvature condition (12.41) the classical equations of motion for the generalized 
NLS model are entailed i.e. 



.dilJi{x,t) d^ipi{x,t) 



%- 



dt d'^x 



+ 2Kj2\i^ji^, t), i, J e{i,...,N- 1}. (3.6) 



It is clear that for N = 2 the equations of motion of the usual NLS model are recovered. 

As already mentioned to obtain the local integrals of motion of NLS one has to expand 
T (12.31) in powers of [22]. Let us consider the following ansatz for T as |A| — oo 

Tix, y, A) = (I + Wix, A)) exp[Z{x, y, A)] (I + Wiy, X))-' (3.7) 

where W is off diagonal matrix i.e. W = Xl j^j ^ijEij, and Z is purely diagonal Z = 



Y,f=iZiiEii. Also 



A" ' ^ A'^ 

n=— 1 n=l 

Inserting the latter expressions (13. 8p in (12. ip one may identify the coefficients W^^'' and Z^^'^ 
(see Appendix B for a detailed analysis). Notice that as iA oo the only non negligible 
contribution from Z^'"^ comes from the Z^|r term, and is given by: 

Z^^lm -L) = iL5n,^r + ^^(^) ^^Ni^)■ (3-9) 

i=l ^-^ 

It is thus sufficient to determine the coefficients wj^^^ in order to extract the relevant local 
integrals of motion (see also [SI])- Indeed solving (12. ip one may easily obtain: 

M/.«(x) = -i^^x), W^Six) = ^^',{x) 

wl^\x) = z^/^i;'^{x)-iJY,\M^)\''Mx)^ •••• (3-10) 
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The Poisson structure for the generahzed NLS model is defined as: 



r \ _ ^ ( 5A SB 6A SB \ 

I ' r '^II U^.(x) <5v^,(x) sux) SM^)^ ^ ' 
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From the latter formulae (13.101) and taking into account (13. 7p . (13. 9p the local integrals of 
motion of NLS may be readily extracted from In trT[L, — L, A), i.e. 

»L N-l 



Ii = —IK / dx 'tpi{x)lpi 

J — hi J 1 



K 

'2 

-IK 



, , .X) 

i=\ 

L N~l 



dx ^ (^i{x)ip[{x) - i)i{x)ip[{x) 

1=1 
L N-l 



I rfx + ... (3.11) 

"^-^ i=l k 



The corresponding familiar quantities for the generalized NLS are given by: 

/i lo U 

u = -—, v = --, n = --, (3.12) 

iK IK IK 

and apparently 

{n, V} = {n, Af} = {AT, V} = 0. (3.13) 

Again in the special case where N = 2 the well known local integrals of motion for the usual 
NLS model on the full line are recovered. 



3.2 The generalized NLS on the interval 

After the review on the NLS on the full line we can come to our main concern, which is the 
evaluation of the integrals of motion after implementing integrable boundary conditions. We 
shall investigate subsequently both SP and SNP boundary conditions. 



3.2.1 SP boundary conditions 

Let us first consider the NLS model with SP boundary conditions. For this purpose c-number 
solutions of the classical reflection equation are needed. A general non-diagonal K matrix 
satisfying the classical reflection equation is given by (see also [6]) 

K{X) = D + A + i^X~^ 

N-l 

D = -Ell -cJ2Eii + Enn, a = 2k{EiN + Em), c = 4k^ + 1. (3.14) 
Apparently in the case where k = a diagonal solution is recovered 

N-l 

K{X) = -Y,En + ENN + i^X~\ (3.15) 

i=l 
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The more general diagonal K matrix is given by (see also e.g. [3l [5]) 

I N 
i=l i=l+l 

The solution fl3.15p may be seen as a special case of fl3.16p for / = — 1. 

Henceforth we set x = 0, y = —L, and we focus on the case with diagonal boundaries 
provided by (IXTSD . (I3T6|) . and K = 1. Also K+{X) = K{X, ^+), ^-(A) = A^(-A, T)- The 
quantity under expansion is 

In tr|A'+(A)T(0,-L,A)K"(A)f (0,-L,A)} = 

In tr|(l + #(0))^i/^+(A)(l + Vt^(0))e^(°-^)(1 + W{-L))-'K-{X){1 + Vt^(-L))e-^(°-^) } 

(3.17) 

where the objects with 'hat' are simply the same as before but now A —A. The technical 
details of the relevant computations are presented in Appendix C. 

(I) Let us first consider the simple boundary conditions described by ( 13.15^ . Gathering 
all the information provided by equations flC.2p . and by explicit computations concerning 
the zA — > cx) expansion (see Appendix C, case (b) for fore details) we conclude that the 
integrals of motion for the NLS on the interval are given by: 

pO N-1 

Ii = —2iK / dx 'ipi{x)ilJiix) , 

J-L -^-^ 

^-1 N-1 

i=l j=l 

N-1 N-1 

+2^^+/€X^/',(0)V',(0)-2ze-«:X^.(-L)V>,(-L), ... (3.18) 

i=l i=l 

the quantity I2 as expected is trivial, as in the case of sine-Gordon model on the half line. 
Recall that in the whole line the quantity I2 corresponds essentially to the momentum, which 
is not a conserved quantity any more. To obtain the number of particles and the Hamiltonian 
we simply have to divide by — 2zk 

A/'=- — , n = -— and {n,Af} = 0. (3.19) 
2iK 2iK 

It is clear that different choices of boundary conditions lead to distinct boundary contribu- 
tions to the integrals of motion. A more detailed description of complicated diagonal and 
non diagonal boundaries is presented in Appendix C. Notice also that for N = 2 the bound- 
ary Hamiltonian presented in [2] is recovered. Of course we could have considered Shcwartz 
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boundary conditions at x = —L i.e. L), ip^—L) = and trivial right boundary K~ = I 
(that is the system is considered on the half line), then the boundary terms appearing at 
X = —L in the expressions of the integrals of motion would disappear. 

As already mentioned the equations of motion will be derived based on the existence of 
a boundary Hamiltonian rather than on the existence of a modified Lax pair. In general, 
among the integrals of motion there exists a Hamiltonian (13.181) such that the relations below 

^^^ = {7^(0,-L), Ux,t)], ^^^ = {7^(0,-L), Ux,t)], 
~L<x<0 (3.20) 

give rise to the classical equations of motion. Indeed considering the Hamiltonian Ti. fl3.18p . 
f l3.19p . we end up with the following set of equations 

z G {1,...,A^-1}. (3.21) 
In general the boundary Hamiltonian for the generalized NLS model may be expressed as 

n= c^o; ^(«:|^,(x)|2^|V',(x)|2 + V':(x)^:(x))+i3 (3.22) 
•^-^ i=i j=i 

where B is the boundary potential. One may write the equations of motion for a generic 
boundary potential B. It is clear that the bulk part remains intact as in fl3.2ip . and what is 
only modified is the boundary conditions at x = 0, x = —L depending naturally on B, i.e. 

d^i{x, t) ^ dB_\ ^ / dij,{x,t) ^ dB_\ ^ ^ ^3) 



dx (9^j/x=o V dx dilJi/x=-L 
Two more examples of diagonal boundaries are presented below: 

(II) Consider the boundary conditions described by f l3.16p . The corresponding contribu- 
tions to the integrals of motion due to the presence of non trivial boundaries are computed 
in Appendix C, case (b). In this case the boundary potential (see flC.6P ) is given by 

N-l 1+ 

^ = - E (v^.(o)^:(o)+v^:(o)^.(o))-e^ 5^ ^.(0)^.(0) 

i=l++l i=l 
N-l I- 

+ (^.(-^)^:(-^)+^:(-^)^^(-^))+r5Z^^(-^)^.(-i^), (3.24) 

i=l-+l i=l 
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and consequently the boundary conditions for the equations of motion at x = 0, x = —L 
now read as: 

^^.+ (0) = ij,-{-L) = 0, j^e{l^ + l,...,N-l}. (3.25) 

The previous case (I) may be seen as a special case of the more general diagonal boundary 
conditions by setting = N — 1. Ultimately, one would like to investigate the SP boundary 
conditions in the context of affine Toda field theories, something that has not been achieved 
up to date. In this case, it is naturally anticipated that the corresponding equations of 
motion should explicitly depend on the parameters contrary to the case analyzed 

in [9], where no extra free parameters associated to the boundaries occur. It is also worth 
stressing that in the context of integrable spin chains the integers appear explicitly in 
the corresponding Hamiltonian as well as in the associated symmetry of the model. More 
precisely, it was shown in [S] that the open spin chain with diagonal boundary conditions 
associated to integers = / is gli ® qIn-i invariant (or Ug{gli) ®Uq{glN^i) invariant in the 
trigonometric case). The symmetry breaking for the quantum gl^ NLS model due to pres- 
ence of non trivial integrable boundaries is also discussed in 



(III) Finally we consider the case where K"^ = I (Appendix C, case (c)). The boundary 
potential in this case is 

^ = - E (^^'(o)^Ko) + Momo)) + E {i^ii-mi-L) + u-Lm-L)) (3.26) 

i=l i=l 

and apparently we end up with simple Dirichlet boundary conditions 

V^i(O) = ^i(-L) = 0, ie{l,...N-l}. (3.27) 

Note that the N = 2 case in particular was investigated classically on the half line in [32] , 
whereas the NLS equation on the interval was studied in 



3.2.2 SNP boundary conditions 

Recall that in this case the object under consideration, compatible with the underlying 
algebra, that is the classical version of the twisted Yangian, is 

In tr|K+ T(0, -L, A) A'" VT\0, -L, -X)V^ (3.28) 

and we choose here for simplicity = I. Note however that a generic solution of the 
classical twisted Yangian is given by any matrix K = (see [H]). We shall choose in 
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what follows V = antid(l, . . . , 1). By expanding (13.281) in powers of A~^, along the lines 
described in Appendix C, explicit expressions for the integrals of motion are entailed (see 
fIC.lSp . flC.16P ). It is worth pointing out, bearing in mind expressions f lC.15|) . that non- 
local contributions to the integrals of motion arise. This is quite an intriguing fact and it 
definitely merits further investigation, which however will be undertaken in a forthcoming 
work. Nevertheless, based on the formulas fIC.lSp . flC.16p we may explicitly express the first 
non-trivial conserved quantity, which is somehow a 'modified' number of particles, i.e. 



/U ru 
dx 'ipi{x)xpi{x) + / dx iJi{x)ipi{x). (3.29) 

Notice that the SNP boundary modify dramatically the number of particles (see (13. lip . 
(I3.12P ). Indeed, the variation due to the integrable boundary conditions is not limited to 
the addition of certain boundary terms to the bulk quantity, as is customary, but it gives 
rise to an alteration of the bulk expression itself. This is a very interesting and definitely 
non-conventional aspect that has not been encountered before, especially in the context of 
continuum integrable theories. Note finally that in the special case N = 2 the 'modified' 
number of particles reduces to the usual number of particles, which is a conserved quantity 
for the s/2 NLS model with diagonal boundary conditions (see (13.180 ). 



3.3 The NLS model with reflecting impurity 

A physically relevant example will be discussed in what follows. More precisely we shall 
restrict our attention to the usual NLS model, and within the framework described in the 
previous section we shall examine the problem of reflecting impurities attached to the ends 
of the system. According to [2] one may consider a more general solution of the reflection 
equation. Consider the classical Lax operator satisfying 

|Li(Ai), L2(A2)} = [ri2(Ai - A2), Li(Ai)L2(A2)], (3.30) 

recall r(A) is given in (13.10 . For example consider the L operator associated to the classical 
Lie algebra s/2: 

LA = ^ 3.31 

' ' \Si+lS2 X-S3 J 

where apparently Sa obey 

{Sa, Sb} = -i'^SabcSc (3.32) 
i=l 
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e being the usual antisymmetric tensor. One may easily express the later matrix in terms of 
canonical variables (x, X), i.e. 

Degenerate cases of the matrix above are for instance the Toda chain and the DST model 
(see e.g. [31] and references therein) with Lax operators given by 

L™"(A) = ( ~^ ) . L--(A) = ( ' ^''^ ) (3.34) 

Consider the following generating function of the integrals of motion 

Intr |K+(A)r(-L,0, A)K-(A)t(-L,0, A)} (3.35) 

K is a generic 'dynamical' type solution of the classical reflection equation [2], i.e. 

K±(A) = L(A - 6) K^{X) h-\-X - 6) (3.36) 

L can be any solution of (13.301) . are any c-number solutions of classical reflection equation. 
Note that the Poisson brackets for K in the classical reflection equation are considered with 
respect to the canonical variables (x, X). Here we shall deal with a simple example, that is 
K"^ = I and 6^ = and L given by (I3.33P (for simplicity set p = 0) then it is clear that 

Finally the boundary contribution to the Hamiltonian is given by (see also Appendix B) 
4''^ = Ihl - hih + /i3 + \hl - hji2 + h + 2iKi){{})i)'{{}) - 2iKij{-L)ij'{-L). (3.38) 

ho = 1, hi = 2Z+, h2 = 2kV'(0)V^(0) - 2iv^(^x+Z+V'(0) + (x+)-^Z+'0(O))) , 

h = 2tKiJj{0yij{0) + 4kZ+^/^(0)^/^(0) + 2v^(^(x+)~^Z+V''(0) + x+Z V(0)) (3.39) 

where Z^ = x^X^. Analogous expressions to fl3.39p are given for h, in particular /i„ = 
(— )'^/i„: —L, {x~^, X+) (x~, X~). Based on the latter expressions the boundary 
part of the Hamiltonian may be deduced. Indeed, bearing in mind that the boundary 
potential is given hj B = — ^ and taking into account fl3.38p . (13.390 we conclude 



i3 = b(x+,Z+,0) -b(x',Z-,-L) (3.40) 
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where we define 

h(x,Z,x) = — '^(xZ'^ijix) +x^^Z'^i!{x)) ^(x~^Zip'(x) +xZ^'(x)) 

y/K\ / i^/K\ / 

- (f{x)ij{x)+tp{x)4j'{x))-—Z^ (3.41) 
V / 3m 

and as expected the boundary contribution of the Hamihonian is solely expressed in terms 

of the canonical variables x^, as well as the boundary values of the fields and their 

derivatives (see also [35] for a similar treatment of the classical sine-Gordon model). 



4 A quantum lattice version of NLS 
4.1 Review on periodic lattice NLS 

Let us first present the general algebraic framework associated to the discrete quantum 
version of the NLS model, introduced and studied for the periodic case in [28l [29]. In the 
quantum level the key object as is well known is the L operator satisfying: 

Ri2{Xi - A2) Li„(Ai) L2„(A2) = L2„(A2) Li„(Ai) R^iXi - X2) (4.1) 
where the i?-matrix associated to the gl^ Yangian is 

R{X) = A - inF, (4.2) 
and obeys of course the Yang-Baxter equation [36l [37], i.e. 

RuiXi - A2) i?i3(Ai) i?23(A2) = i?23(A2) R^iXi) Rui^i - A2). (4.3) 

We shall focus in this and the subsequent section in the simplest s/2 Yangian. In this case 
a simple solution of equation (14.11) is given by (on a detailed description of the underlying 
algebra see e.g. [281 [29] ) 

^ M^ _ I 1 - «AA + A2k0„^„ -iA^(pn , 



where the ipn, (pn satisfy canonical commutation relations 

i'n, <t>, 



^Snm- (4.5) 



In fact this solution may be thought of as the quantum version of the NLS. Indeed the 
classical limit of the L operator (14. 4p gives U ( 13. 3p (for further details see [28]). Set V'n = 



j^„+A ^(^x) then as A ^ 

^fA^V where fK.T. = I 

^/kiI){x) 



L(A) = l + AU(x,A) + 0(A^), where U(x, A) = I ^ \f^(t>{x) . ^^^^ 
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note that (f){x) = ipix), and U is equivalent to U (13. 3p of NLS up to a gauge transformation 
i.e. 



(4.7) 



It is more convenient for our purposes here to use L with a rescaled spectral parameter 
matrix. Let us multiply (14.41) by ^ and also set C = f then the rescaled L matrix may be 
written as: 



Un(A) - I ^ 

N„ = 1 + K^^cpni^n- For the special value C = the L operator reduces to a projector 



(4.8) 



L(C = 0) 



1 




(4.9) 



Due to the fact that the algebra (14. ip is equipped with a coproduct one may build tensorial 
representations and construct a spin chain like system with periodic boundary conditions, 
by introducing the quantities 

To(A) = LoL...Loi, and t{X) = tr^ Tq{X) (4.10) 

with T{\) being essentially the quantum analogous of (12.31) and apparently by virtue of (14.11) 



t(A), t(A') 



0. 



(4.11) 



By expanding Int(C) around C = we find the corresponding involutive quantities exactly 
as in the classical case |28]. It is easy to see that due to (14. 9 p t(0) = 1, then one finds (for 
more details we refer the reader to |28j ) 



Ini(C) = C^C*! + C'kCs + C'/^C-s + 







(4.12) 



with 



n=l 

L L ^ 

n=l n=l 
^ L L 

^^ = ^H^- = 11 [^n+lll^n-l - (N, + N,+l)0n+lV'n + (3kA2)-1N^) . 



(4.13) 
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From the latter objects one may derive lattice versions of the classical quantities (13.121) . 

TV = (Ci - L)|a^o ^ dx (j){x)ip{x), 



V = {C2 + ^^)Ia^o dx {<l>.,i) - (P^,) 

L 



n = -C; + ^;^\a^o ^ J dx (0,^, + k(0^A)2). (4.14) 

Notice that the expressions above are symmetric to "0, so one can set = -0 and obtain the 
familiar expressions for the NLS system (13.121) . Note also that the existence of an obvious 
pseudo- vacuum allows the implementation of Bethe ansatz techniques [38| [28]. however our 
aim here is to extend such computations in the case of the integrable open spin chain, which 
is discussed in the subsequent sections. 

4.2 Open lattice NLS 

We come now to the case with open boundary conditions, which is our main concern. The 
underlying algebra in this case is defined by the reflection equation pQ 

RuiXi - A2)/^i(Ai)i?2i(Ai + A2)/^2(A2) = K2(A2)/22i(Ai + A2)i^i(Ai)/2i2(Ai - A2). (4.15) 

The tensorial type solutions of the reflection equation as well known is given by [2] 

%{X)=To{\) K^{\,r,c') T,-\-X) (4.16) 

K~{X,^~ ,k~) is c-number solution of the reflection equation with the most general form 
given by ^9] 

i^±(A) = Xa' ± i^^ + 2A;^A(a+ + a'). (4.17) 
Note that the explicit expression of L~^(— A) = ]L(A) is given by: 

lon{X)=( ' ^^^'^^ V (4.18) 

The corresponding generating function of the conserved quantities of the open system is 

tiX)=tro{K^iX,e,k^) ToiX)} (4.19) 

K~^{X, k~^) = K{—X—i, , k~^y, K is again a c number solution of the reflection equation. 
And due to (I4.15P it is clear that integrability of the quantum system is ensured i.e. 

t(A), t(A')l = 0, A, A' G C. (4.20) 
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In the remaining of this section we specify the spectrum of the lattice NLS model with 
diagonal boundary conditions by means of the Bethe ansatz technique [38]. Focusing on 
diagonal boundaries should be sufficient given that in ^ the spectral equivalence between 
systems with diagonal and non diagonal boundaries was shown by means of appropriate 
gauge transformations, but only for spin chains associated to the fundamental representa- 
tion of sl2- Presumably there exist suitable gauge transformations for the system under 
consideration, such that the spectral equivalence is guaranteed. We shall further comment 
on this point on a separate publication. When both boundaries are diagonal there exists an 
obvious reference state for the transfer state 

N 

fl = ^'cun with ^„ ro„ = 0. (4.21) 

n=l 

Based on this observation one may in a straightforward manner derive the spectrum and the 
corresponding Bethe ansatz equations. We provide directly the results avoiding the technical 
details (for a detailed description we refer the reader to [2]). The spectrum of the transfer 
matrix is given by 

A(A) = g(A)b,(A) n + + h(A)b.(A) n (A-A.-y(A + A,-2..) 

(A-A,)(A + A,-z«:) ^ (A-A,)(A + A,-z«:) 

(4.22) 

where we define 



g(A) = i-iXA + 1)^, h(A) = (zAA + 1 + kA)^, 

bi(A) = |^(A + ^^)(-A + te), b2(A) = -^(-A + + ^K){X -tK + z^) 

(4.23) 

The corresponding Bethe ansatz equations arising as analyticity conditions on the spectrum 
are: 

g(A» + f) bi(A, + f) _ TT Ai - Aj - in Xj + Xj - in (AOA) 
HXi + f ) b2(A, + f ) ~ f}^ Xi - X, Xi + Xj + tK ^ ' ' 

Notice that although we deal with an open spin chain and one would expect a leading order of 
2L, we see a leading order of L exactly as in the periodic case. The same phenomenon occurs 
in the boundary lattice Liouville model [ID], and is presumably associated to the degenerate 
nature of the L matrix (see also similar comments in [M]). It should be emphasized that the 
Bethe ansatz equations (14.241) are of particular significance given that their thermodynamic 
analysis yields for instance consequential information regarding the corresponding bulk as 
well as boundary exact S matrices of the model. 
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5 The generalized lattice NLS 



We shall deal in what follows with the lattice quantum version of the qIn NLS model. Recall 
that the gl^ Yangian R matrix, solution of the Yang-Baxter equation (14. 3p . is given in (14. 2p . 
The relevant L operator in this case is given by (see also |28j ) 

N~l N N 

L(A) = {-- + Y,<P^H^^^)E,, + Y,E,, + Y.{<P^^~'^E,, + ij^^~'^E^ (5.1) 
i=i i=2 j=2 

Notice that here we set implicitly iA^/K = 1. It will be also useful for the following to define 

UX) = ViV'i~\ + tKp)Vi, p = j (5.2) 
we choose here V = antid(l, . . . , 1), which gives rise to the following explicit form: 

N-1 N-1 N 

L(A) = C- + p+J2 <P'^H^'^)Emn + E,, + Y,{<P^'-'^E-,^ + ^^-^-^^E^^) (5.3) 
j=i j=i i=2 

where j = N — j + 1. Recall that in general we could have chosen any V such that V'^ = I. 
We shall focus hereafter in the case of SNP boundary conditions, given that they are not so 
widely known compared to the SP ones, especially in the context of integrable lattice models. 
The main objective in this section is to derive the exact spectrum and the corresponding 
Bethe ansatz equations. Note that in the SNP case the underlying algebra is defined by the 
following relation, (twisted Yangian, see e.g. [T7] ) 

i?i2(Ai - A2) K^{X^) i?2i(Ai + A2) i^2(A2) = i^2(A2) RuiXi + A2) A'i(Ai) i?2i(Ai - A2)(5.4) 
and in analogy to the classical case we define 

R{X) = Vi R{\{-X + ipn) Vi. (5.5) 
The generating function of the quantum integrals of motion in this case is defined as: 
t{X) = tr |i^+(A) T(A) A'~(A) T(A)}, with 

T(A) = Lol(A) . . .Loi(A), T(A) = Loi(A) . . .Lol(A) (5.6) 

are c-number solutions of the twisted Yangian (15. 4p . In fact, it was shown in [13] that 
any matrix K = is a solution of the twisted Yangian. In Appendix D an explicit 

computation of local integrals of motion for particular boundary conditions is presented. 
Based on these findings we present the explicit form of the boundary momentum in the case 
where = V, i.e. 

L L-1 N-1 N-1 

= -y ( E - 2 E E ^^^^ + + • (5.7) 

n=l n=l j=l j=l 
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We could have of course considered K"^ oc I, which is the case will be examined subsequently, 
but for the sake of simplicity we considered the aforementioned boundary conditions, which 
from a technical point of view are much easier to deal with. By comparing the bulk mo- 
mentum given in Appendix D flD.4|) with (15.71) we conclude that the periodic terms in flD.4l) 



are replaced essentially by the last two boundary terms in flD.7l) . whereas the bulk part re- 
mains intact. Following the logic of fl4.14p . it is expected that the expression (15.71) should be 
a regularization of the continuum generalized NLS model momentum with particular SNP 
boundary conditions, exactly as it happens for the periodic NLS model (indeed compare the 
bulk continuum expression (13.111) with the discrete periodic analogous (1D.6P ). Comments 
on higher conserved charges may be also found in Appendix D. 

To deduce the spectrum and Bethe ansatz equations of the generalized NLS model with 
SNP boundary conditions we shall restrict our attention to another simple case i.e. K"^ = I. 
Again the spectral equivalence between systems with diagonal and non-diagonal boundaries 
is discussed in [HI [16], for spin chains associated to the fundamental representation of gli\[. 
The first step toward the diagonalization of the transfer matrix ( 15.60 is the derivation of a 
reference state. Indeed, in this case there exists an obvious reference state, that is 



L 



n = 6dwn: iJ^^Wn = 0, ne{l,...,L}, t e {1, . . . N - 1}. {51 



n=l 



The corresponding eigenvalue may be easily derived using the fact that L, L and conse- 
quently T, T satisfy 

Ti(Ai) i?i2(Ai + A2) T2(A2) = T2(A2) i?i2(Ai + A2) Ti(Ai). (5.9) 

Taking into account the latter relation we conclude that the actions of the transfer matrix 
on the pseudo vacuum provides the following eigenvalue: 

N-l 

A(o)(A) = a^(A) g^{X) + ^ <^„(A) + F(A) g^{X) (5.10) 

n=2 

where we define 

/s\ ^A iX 
a{X) = , b{X) = — + p 

K K 

A-f(p-l) N + 1 

9niX) = ^^-j^ , 1 < n < 

Q iv+i =1, odd 

A + 1 

gi = gTi-X + iKp), ——<l<N. (5.11) 
The functions gn are essentially 'boundary' contributions to the spectrum. 
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To determine the general eigenvalue form we shall adopt the analytical Bethe ansatz 
formulation [H]. The basic assumption within this framework is that the structure of any 
eigenvalue is similar to the pseudo-vacuum eigenvalue i.e. 



AT-l 



A(A) = a^(A) ^i(A) A(A) + J] g^{\) A^{X) + h\\) (?^(A) A^(A). (5.12) 



n=2 



The so called dressing functions An may be explicitly determined by imposing certain phys- 
ical and algebraic requirements, such as analyticity, crossing, etc. We do not give all the 
details of the formulation, (for a more detailed description of the process we refer the reader 
to [U [6l [H]), but the explicit expressions for the dressing functions are given by: 

A + \f + f A - \f + f 



^i(^) = n . , ,(1) . , ,(1) 



A + A • - Y A - A . - Y 



Y -ti^ A - A^'^ - 2 



X + X^P X- XP -^^-iK 



\ I \(0 i^K \ \(0 Uk 

M('+^)^^um)_^^« A-A^^^-^ + f . N-l 



n l—l 2 ' 2 2 a 2 ' 2 l</< 
\ , ilK iK \ Uk ift ' — 

A + A^. - — - y ^ - - — - T 

- 1 

Ai{X) = Ai{-\ + iKp), -——<1<N (5.13) 
and in particular for = 2n + 1 

An+l(A) = TT ^ ^ 2 

7 J_J_ \ I \C«) iriK \ \{n) iuK 

i=l ^ + ~ ~ ^ ~ ~ ~ 

M("+1) ^ I iriK I IK \ iUK , ik 

11 X I _ inK _ iK \ _ _ voK _ ik' \ ' J 

j=l ^ + 2 2 ^ 2 2 

Finally Bethe ansatz equations follow as analyticity requirements upon the spectrum, and 
they are written explicitly as: 

(i) N=2n+1: 

("(a!" + f )) in + (1 - Sa) = - n <^M' - )e.(A« + AJ) 
X n -,(Af-Af")e_,(Af + Ar') 

/ = 1 77,-1 
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i(Af ^) = - n ^2(a!") - ))e2(Af ) + )e„i(Af ) - Af ))e„i(A;") + Xf^) 



X 



n e_i(A;") - Af -^))e_i(AS") + X^;-'^) 



(5.15) 



with r 



±1, M(^+i) = and define e„(A) 




(ii) N=2n: In this case the Bethe ansatz equations for / = 1, . . . n — 1 are the same as 
in the previous case. What is only modified is the last set of equations, which takes the 
form: 



e-i(Af ^) = - n ^^(^' - Af ^)e2(AS") + Af n ^-i(Af^ - Af -^V^_,(a!") + Af-^Y 



Such type of boundary conditions were first introduced in [13] for the sis spin chain, whereas 
generalizations investigated in [6],[T3] from the physical (Bethe ansatz) as well as the algebraic 
point of view. The associated symmetries were studied in detail in [HI [151 IIS]) while in 
[T5] both SP and SNP boundary conditions were examined in parallel. The interesting 
observation is that the RHS of the equations above in the case where = 2n+l coincide with 
the ones associated to the osp{l\2n) algebra. In any case the Bethe ansatz equations (15.151) 
are somehow 'folded' compared to the usual gl^ ones. This is expected given that folding 
occurs at the algebraic level as well (Dynkin diagrams), and only the subalgebra invariant 
under charge conjugation survives after the implementation of these rather unconventional 
boundary conditions (see also relevant comments in [H]). The case of SP boundaries can be 
also treated along the same lines, and the entailed spectrum and Bethe ansatz equations will 
have the usual gl^ structure (the expressions are omitted here for brevity). More precisely 
the LHS of the Bethe ansatz equations will have exactly the same form as the usual gl^ BAE 
for an open chain (see e.g. [311S11S]), while the RHS will depend explicitly on the actions of 
the diagonal entries of the L, on the pseudo-vacuum. 

6 Discussion 

To summarize, SP and SNP boundary conditions were studied for the classical generalized 
NLS model, and the boundary integrals of motion as well as the relevant classical equations 
of motion were explicitly derived. This was the first time, to our knowledge, that SNP 
boundaries were implemented in the context of the generalized NLS model. Nevertheless, 



(5.16) 
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there are still several open questions especially regarding the locality of some of the integrals 
of motion for particular choices of left/right boundaries, which however will be left for 
future investigations. In the same spirit the usual (s/2) NLS model with reflecting impurities 
was also analyzed. Moreover, a suitable lattice version of the generalized NLS model was 
considered and the SNP boundary conditions were implemented, given that in general they 
are much less studied in this context. For this choice of boundary conditions we were able to 
specify the first non-trivial local integral of motion i.e. the 'boundary momentum'. We also 
derived the spectrum and Bethe Ansatz equations for the simplest left/right boundaries. 

Although SP boundary conditions are somehow the obvious ones in the framework of lat- 
tice integrable models, they have not been really considered up to now in classical continuum 
integrable theories. Therefore, it will be our next goal to impose the SP boundary conditions 
to other well known classical systems such as (massless) affine Toda field theories, principal 
chiral models, etc. In addition, the investigation of the generalized (m)KdV hierarchies [12] 
with integrable boundaries is another very interesting direction to pursue together with their 
quantization into an appropriate lattice version (see e.g. [43] and references therein). Once 
this point is clarified, the study of the underlying dynamical symmetries constrained by in- 
tegrable boundary conditions could be discussed in full generality along the lines described 
in [H] , and this would definitively shed new light on the character of the different integrable 
boundary conditions. More precisely, it would be of great consequence to examine how the 
so called hidden symmetries constructed in [H] are modified in the presence of non-trivial 
integrable boundaries, and in particular in the case of (quantum) twisted Yangians. Finally, 
an interesting direction to pursue is the explicit derivation of the modified Lax pairs by 
means of the 'boundary' integrals of motion. We hope to report on all these issues in forth- 
coming publications. 
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number HPRN-CT-2002-00325. 

A Appendix 

In this appendix the classical integrability for models with both types of boundary conditions, 
SP and SNP, is reviewed. The first step in order to prove the classical integrability is to show 
that the quantities introduced in (12.111) are indeed representations of the algebras defined 
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by (12.91) . To achieve this we shall need in addition to (12.61) the following set of algebraic 
relations emerging essentially from (12.61) . i.e. 

{ri(Ai), f2(A2)} = ri2(Ai - A2) Ti(Ai) r2(A2) - Ti(Ai) T2(A2) ri2(Ai - A2) 

{ri(Ai), f2(A2)} = Ti(Ai) r^X, + A2) T2(A2) - T2(A2) rl,{X, + A2) Ti(Ai) 

{ri(Ai), T2(A2)} = Ti(Ai) rUX^ + A2) T2(A2) - T2(A2) r*2(Ai + A2) Ti(Ai) (A.l) 

Our aim now is to show that (12.91) are satisfied by (12. 8p : 

{Ti(Ai), T2(A2)} = [T^{X,)K^{X^)f^{X^), T2(A2)i^2"(A2)T2(A2)} = ... 

= Ti(Ai)T2(A2) ( - Ki(Ai)K2(A2)ri2(Ai - A2) + ri2(Ai - X2)K^{X^)K2iX2) 

+Ki{Xi)rUXi + X2)K2{X2) - K2(A2)rt2(A2)i^i(Ai))ri(Ai)T2(A2) 
+Ti(Ai)r*2(Ai + A2)T2(A2) - T2(A2)r*2(Ai + A2)Ti(Ai) 

ri2(Ai - A2)Ti(Ai)r2(A2) - ri(Ai)T2(A2)ri2(Ai - A2) (A.2) 

and making use of (12.101) . (1A.3I) and ( lA.ll) we end up to (12.91) . Recall also that c-number 
solutions of the above equations satisfy the following 

ri2(Ai - A2), J^i(Ai)K2(A2)] = K2{X2)rUX^ + A2)i^i(Ai) - K,{X^)rl,{X^ + A2)/^2(A2), 

(A.3) 

which is equivalent to (I2.10p . 

We may now show exploiting (12.91) and (1A.3P the classical integrability (12.131) . Indeed 
consider the following object 

{aY(Ai)Ti(Ai), /^2^(A2)r2(A2)} (A.4) 

then taking the trace in both spaces 1 and 2, and considering the defining relations (12. 9p . 
(I2.10p we end up with 

|t(Ai), t(A2)} = tri2(Ki+(Ai)K2+(A2)ri2(Ai- A2)Ti(Ai)r2(A2) 

- i^+(Ai)A'2+(A2)Ti(Ai)T2(A2)ri2(Ai - A2) 
+ K+(Ai)ir2+(A2)ri(Ai)r*2(Ai + A2)T2(A2) 

- K+(Ai)At(A2)r2(A2)rt2(Ai + A2)Ti(Ai)). (A.5) 



Finally moving appropriately the factors of the products within the trace and using (]A.3P it 
is straightforward to show 

|t(Ai), t(A2)} =0, Ai, A2eC (A.6) 

and this concludes our proof. Similar arguments hold also in the case of dynamical boundaries 
discussed in section 3.3. 
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B Appendix 



We present here some technical details on the derivation of the conserved quantities for 
the generalized NLS model on the full line. The first step is to insert the ansatz (13. 7p in 
equation (12. ip . Then we separate the diagonal and off diagonal part and obtain the following 
expressions: 

Z' = AUi + {VqWY^^ 

W + WZ' = l]o + {VoWY^^ + XUiW (B.l) 

where the superscripts (D), (O) denote the diagonal and off diagonal part of the product 
UoM^. Recall that W = ^ijEij, Z = J2i ^uEu then it is straightforward to obtain: 

7V-1 
i=l 



i=l 

{U,W)^o) J2 [iPWnjE,^ + 'l^J,W^JEN^^ . (B.2) 

iy^j, iy^N, j^N 

Substituting the latter expressions flB.2l) in fIB.ip . we obtain 

N-l N~l „L 

Z{L, -L, A) = -tXL(^Yl - ^^^) + / dx(ij,WmEn + ^iW,^£^jv7v) (B.3) 

1=1 i=i 

And recalling that the leading contribution in the expansion of (In trT) , -T is given in (13. 7p - 
as iA — > oo is coming from the ^atat term we conclude: 

Znn{L, -L, A) = iXL + VK.y^l dx 'ilJi{x)WiNix) (B.4) 

1=1 ^-^ 



Due to ( ]B.4p it is obvious that in this case it is sufficient to derive the coefficients WiN only. 
In any case one can show that the coefficients Wij satisfy the following equations: 



W[jE,, -iXY^ (^NiEm - W^nE^n) + v^J^ (^JiWl.Em + ^^W^J,E,M) = 
Z (^i EiN + ipiEm^ + ^/K Y (^iJi Wnj Eij + tpi Wij E^j ) 
-v^ J2 (i^jWNjWi.Eij+^iWiNWjNEjN) (B.5) 



iy^j, if^N, j^N 



Finally setting Wij = Xl^i using (IB.SP we find expressions for W-j^ reported in 

(I3.10p . In the case with integrable boundary conditions, we shall need in addition to (13.101) 
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the following objects: 



iy^j, j^N 



W[f = tWll^^W;;] - ^W^(1W]<^ (B.6) 

C Appendix 

In what follows we evaluate the boundary terms contributing to the Hamiltonian for right 
and left boundary described by the more general, diagonal and non-diagonal, solutions of the 
reflection equation (SP boundary conditions). Moreover, for the SNP boundary conditions 
we identify the corresponding integrals of motions, and we explicitly evaluate the first non- 
trivial charge. 



C.l SP boundary conditions 

We shall expand the generic object (13.171) keeping of course only diagonal contributions. 
More precisely, as in the bulk case due to the fact that the leading order is e*'^'^ as — oo 
the only non negligible part is coming from the Enn terms, hence we shall only consider 
such contributions: 

(1 + 1^(0, X)y'K+i\)il + WiO, A)) 

n=0 

"(1 + W{-L, X))-'K-{X){1 + W{-L, A))j = ^ 

L J NN ^ — ^ A 

n=0 

Z(0,-L,A) - Z{0,-L,X) 

n=l 

Again considering only the contribution of the term e^^^ as iA — oo we end up with the 
following expression 

In tr|K+(A)r(0,-L, A)A'-(A)r(0,-L,A)| = 

n=l n=0 n,m=0 




iXL + J2{ 



A" 



(C.l) 
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Recall from (lC.2p that the boundary contribution lies basically in the logarithmic function, 
hence one has to expand the log i.e. 

°° h 4-h °° h h °° f 

n=0 n,m=0 n=0 

where /„, provide essentially the boundary contribution to the integrals of motion (plus 
possible total derivatives from the bulk part) for the left right boundary respectively. The 
interesting observation is that the boundary contribution decouples nicely to terms associated 
to left and right boundary separately, i.e. no mixing occurs 

fi = hi + hi, f2 = -^/ii + /i2 - ^hl + h2, 

h = - + h + Hi - hji2 + h, ... (C.4) 



(a) We first consider generic non diagonal boundary conditions described by f l3.14p . One 
can explicitly evaluate /i„ for the generic case: 

ho = 1, hi = -2ik+y^(^iJi{0) - ^i(O)) + it, 

N-l 

h = 2A;+v^(vi(0) - V'i(O)) - ^k+\ ^ i^^{0)i>^{0) , 

1=2 

N-l N-l 

N-l 

+2ik+J lV^j(0)l'(V'i(0) - 2V^i(0)) + 2tk+K^iijl - 4iKk+^ (^i(O)^j(O) + ^i(0)V'j(0) 

j i-2 



(C.5) 



In fact it is clear from the expansion of the left and right boundary contribution that /i„ = 
: -L, ^+ r, k+ k-. 

Given that there is an overall derivative from the bulk part of Z^, giving rise to a boundary 
term {2iK,Yli'^i'^i), conclude that the boundary contribution to the conserved quantity 
J3 (i.e. the Hamiltonian) is given by 

N-l N-l 

h = h + 2iK Y - 2^«^ Y i''ii-L)M~L) = -2tKB, (C.6) 

i=l i=l 

where B is the boundary potential, and recall that is given by flC.4l) . For diagonal bound- 
ary, terms proportional to k^ apparently disappear, and the Hamiltonian reduces to fl3.18p . 
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For a purely antidiagonal boundary only terms proportional to survive, all the other 
terms disappear. 



(b) The more general diagonal boundary conditions are described by fl3.16p . We associate 
the integers and the free parameters to the right/left boundaries. In this case from 
the expansions fIC.ip and taking into account fl3.10l) . flB.6p we find: 

N-l 

i=i++i 

1+ N~l N-l 

/i3 = -2z«:J]^,(0)V^K0) + 2?«: ^ ^^^0)^,(0) + 2^/^e ^ ^i(0)^/;,(0), ... (C.7) 

1=1 i=l++l i=l 

similar expressions of course hold for hn, i.e. /i„, = (— l)"/i„: —L, — > l~ , 
Taking into account (lC.4p . (1C.6P and derivative contribution from the bulk {^inY^fj^^ 
we conclude that the boundary contribution to the Hamiltonian is given by: 

N-l N-l 

h = 2t^[ (^,(o)^:(o)+^:(o)^.(o))- Y {M-Lm-L)+^:{-m{-L))' 

i=l++l i=l-+l 
1+ l- 

+ 2zK[e+^^i(0)^A,(0)-r$^^^(-^)^^(-^)]- (C.8) 

i=l i=l 

Note that when = — 1 one recovers the diagonal limit of the more general case (a). 

(c) Finally the case where = I may be treated in the same spirit. In this case we 
find that 

N-l N-l 

ho = l, hi = 0, /i2 = 2/€^ ^,(0)^,(0), /i3 = 2z«: J]^,'(0)^i(0), ... (C.9) 

i=l i=l 

and the corresponding boundary contribution to the Hamiltonian is given by 

N-l N-l 

h 



[ Y (^.'(0)^.(0) + ^.(o)^:(o)) - Y (i^^i-mi-L) + u-mi-L))] .(c.io) 



1=1 i=i 



C.2 SNP boundary conditions 

Recall that in this case the object under consideration is given in (13.281) . also we consider for 
simplicity = I and we choose V = antid(l, . . . , 1). Before we proceed with the evaluation 
of the integrals of motion let us first introduce some useful notation 
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°o r(n) 

and (1 + W^(A))-^ = 1 + F(A) where F(A) = V — , (C.ll) 

Z — / \n 

n=l 

recall j = N — j + 1, and also one may easily compute that 

f(i) = -iyW, f{2) = (v[/(i))2_v[/{2), ... (C.12) 

We shall need the following contributions in order to evaluate the corresponding integrals of 
motion: 

(1 + W{0, A))(l + W{0, A)) = 1 + H{X), 

(1 + W{-L, X)y\l + Wi-L, A))-i = 1 + ^(A). (C.13) 

Also bearing in mind that the leading contributions in the considered expansion, as iA ^ oo, 
are coming from Z^vat and for z G {2, . . . , A^} we may write: 

~ Zi;l(0, -L) + (-l)"Zi"Vo, -L) 
ZMO, -L, A) + Zu{0, -L, A) =tXL + J2 \n ' (^-14) 

n=l 

Gathering all the information given above we end up with the following expression 

- 4t(o,-^) + (-ir^fiHo-L) 



In 



tr|r(0, -L, X)VT\0, -L, -A)f} =iXL + J2 



n=l 



N 

+ ln (l + //^^(A) + ^^^(A) + ^e^-(°-^'^)-^--(°-^'^)i/,;v(A)i/^,(A)) (C.15) 

Finally taking into account the information provided above (lC.lip - (lC.15P we may express 
the first non-trivial integrals of motion as: 

h = Z(;^(0,-L)-Z«(0,-L) 

(C.16) 



U l>{J 

dx ■ipi{x)ipi{x) — IK dx ^pl{x)'^jJl{x) 



Notice that in general due to the presence of in ( 1C.15|) non-local terms seem to arise 
in the higher integrals of motion, which is quite an unusual issue and shall be addressed 
elsewhere. Nevertheless, a straightforward computation of the higher charges, based on the 
explicit expression (IC.lSp . may prove the locality or not of the higher integrals of motions. 



Moreover, the presence of Z^^^ alters the structure of the bulk part of the integrals as well. 
The latter integral of motion (1C.16I) gives rise to a 'modified' number of particles. Mm = 
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D Appendix 



In this appendix we shall evaluate the first integrals of motion of the quantum discrete glN 
NLS model with SNP boundary conditions. This model may be also regarded as a higher 
rank algebraic extension of the SI2 DST model (see e.g. [31]), holding a special place between 
the glisf quantum spin chains -extensions of the Heisenberg model- and the qIn generalization 
of the Toda chain. To explicitly specify the local integrals of motions of the model with open 
boundary conditions we shall, as usual, consider the asymptotic expansion of the generating 
function T(A). We shall focus here on the simple case where both left and right boundaries 
are given by K^{X) = antid(l, . . . , 1), and effectively we shall expand 

t{X) = tr T{X) r(A) where f (A) = r*(-A) (D.l) 

and recall T(A) is given by (15. 6p . Indeed after expanding in powers of A~^ we obtain 

T(A) oc En + \t^'^ + ^^^'^ + ^(^"') 
A A^ 

f (A) oc + if + ^f^'^ + OiX-^), (D.2) 
A A 

where the quantities T^^''^\ tO->'^) are defined below 

L N N 

n=l j=2 j=2 

L N N 

n=l j=2 j=2 

L-1 N N 

n>m n=l j=2 j=2 

L-1 TV N L N N 

+ E E ^t'^E,, + Y <^r'' E ^nE,, + y: ^tpE,, + Y <i>t''E,,) 

n=l j=2 j=2 n=2 j=2 j=2 

L-1 N N 

t(^) = -.^[Y ^-^-Eu +Yi: ^r^vti'^i^n + Y ^'"vr^^i?,, 

n<m n=l j=2 j=2 

L-1 N N L N N 

+ E E ^t'^E,, + Y 4>'r' Y ^nE,i + Y ^tpE,, + Y <^t''E,,) 

n=l j=2 j=2 n=2 j=2 j=2 

N-1 

where Nn = Y ^n'^'t'n\ N„ = N„ + p. (D.3) 

In the expressions above all the lower indices denote the site of the spin chain, while the 
upper indices denote the component of the (A^ — 1) dimensional vector fields. 
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We may easily obtain first the bulk integrals of motion by considering the expansion 
tr T(A) = E„ i^- Indeed for the bulk case after simply taking the trace of T^^'^^ we obtain 

L L-l N-l N-l 

n=l n<m n=l j=l j=l 

The quantities identified as the number of particles and the momentum in the NLS model 
are given by the following expressions 

Afd=-h, V, = -(ll!-h) (D.5) 

and more precisely 

L . L L-l N-l N-l 

n=l n=l n=l j=l i=l 

We come now to the open NLS model, and we consider the expansion of fID.ip . The first 
charge of the open model is zero, that is the number of particles is not a conserved quantity 
anymore. The second charge is given by 

L L-l N-l N-l 

n=l n=l j=l j=l 

and corresponds to the momentum 'Pd = which is obviously modified due to the presence 
of the open boundaries. The third charge again is trivial, involving only boundary terms. 
We do not compute any higher conserved charges, but we may rather safely conjecture that 
the only non-trivial conserved charges are the even ones. 
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